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Abstract 

We investigate the interplay between the electron-electron and the electron-phonon 
interaction in the Hubbard-Holstein model. We implement the flow-equation method 
to investigate within this model the effect of correlation on the electron-phonon effec- 
tive coupling and, conversely, the effect of phonons in the effective electron-electron 
interaction. Using this technique we obtain analytical momentum-dependent expres- 
sions for the effective couplings and we study their behavior for different physical 
regimes. In agreement with other works on this subject, we find that the electron- 
electron attraction mediated by phonons in the presence of Hubbard repulsion is 
peaked at low transferred momenta. The role of the characteristic energies involved 
is also analyzed. 
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The electron-phonon (e-ph) interaction can play an important role in many 
physical systems, where it can be responsible for instabilities of the metallic 
state leading, e.g., to charge density waves and superconductivity. When the 
e-ph is particularly strong polaronic effects can also be responsible for a huge 
increase in the effective electronic masses. On the other hand the issue of strong 
electron-electron (e-e) correlations is of relevance for many solid state systems 
like, e.g., (doped) Mott insulators, heavy fermions, and high-temperature su- 
perconductors. Therefore the interplay between the e-e correlation and the 
e-ph coupling is far from being academic and its understanding is of obvious 
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pertinence in several physically interesting cases. In particular, in the super- 
conducting cuprates the e-e correlations are strong and are customarily related 
to the marked anomalies of the metallic state, including the linear tempera- 
ture dependence of the resistivity((l|), (@), where phonons are not apparent. 
This led to believe that the e-ph interaction plays a secondary role in these 
materials, and that the main physics is ruled by electron correlation. How- 
ever, the presence of polaronic spectroscopic features (0) as well as various 
isotopic effects (4;) in underdoped superconducting cuprates shows that the 
e-ph coupling is sizable in these systems. Recent results obtained with Angle 
Resolved Photo-Emission Spectroscopy (ARPES) (jjj have also been inter- 
preted in terms of electrons substantially coupled to phononic modes. These 
contradictory evidences naturally raise the issue of why correlations can re- 
duce the e-ph couplingin some cases and not in other. In this regard it has 
been often advocated (0; H; 0; Hll ) that strong electron correlation suppresses 
the e-ph large angles scattering, which mainly enters in transport properties, 
much more than the forward scattering, which is instead relevant in deter- 
mining charge instabilities (0). Therefore the momentum dependence of the 
e-ph interaction can give rise to important effects and select relevant screening 
processes due to correlations. 

In this paper we study the effect of (a weak) electronic correlation on a two 
dimensional Holstein e-ph system using the flow-equation technique developed 
by F. Wegner ^3 . Ilif ). The flow-equations describe the evolution of the pa- 
rameters of the Hamiltonian under an infinite series of infinitesimal canonical 
transformations labeled by a parameter / that rules the flow. The generator 
of the canonical transformation is chosen in such a way to eliminate some 
interactions in favor of other interactions. Shall one solve the flow equations 
exactly the initial hamiltonian at I = and the final hamiltonian at / = oo 
are guaranteed to have the same eigenvalues. Of course in general the flow- 
equations are not solvable, but suitable approximations can be introduced to 
obtain qualitative trends. Under certain approximations one can map the in- 
teracting problem into a noninteracting one. In this case the method can be 
roughly seen as a way to derive a Fermi liquid theory where bare fermions 
evolve into quasiparticles under the flow. Since one can choose what parts 
of the hamiltonian remain, the method can also be exploited to get effective 
interactions. Of course, in getting a fully diagonal Hamiltonian one is nat- 
urally led to discover the "low-energy" phsysics of the problem. This is the 
approach (and the spirit of it) previously adopted by Wegner and coworkers 
to analyze the instabilities of various forms of the Hubbard model (|l4t ll5tll6|) . 
Another application of the method is the derivation of effective interactions. 
For example this scheme was used to obtain an effective e-e interaction by 
eliminating phononic degrees of freedom in the absence of the initial e-e in- 



teraction fjlTr) . Rather than obtaining the low-energy physics one obtains an 
interacting Hamiltonian which still has to be solved. Analyzing the structure 
of the resulting effective interaction we will discuss the physics of the e-ph 
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interaction in the presence of the e-e interaction and vice versa. We will see 
that this procedure provides the correct energy scales for the phonons dressing 
the effective e-e interaction and, conversely, of the e-e interaction dressing the 
e-ph interaction. 

The paper is organized as follow: first, in Sec. 1, we introduce the Hubbard- 
Holstein Hamiltonian and, in Sec. 2, we present the formalism of the flow 
equations technique. Then in Sec. 3.1 we validate the method by applying it 
to a two atom system. In Sec. 3.2 we compute the effective e-ph interaction in 
a two-dimensional system and in Sec. 4 we eliminate the e-ph interaction to 
get an effective e-e coupling, thus providing a model Hamiltonian suitable to 
discuss different physical regimes. Finally the summary and the conclusions 
are provided in Sec. (5). 



1 Model 



We consider a system of interacting electrons and phonons described by the 
Hubbard-Holstein model, where the electron density is locally coupled with a 
dispersionless phonon. For the sake of simplicity and since this case is relevant 
in physical systems of broad interest like the cuprates and the manganites, we 
consider the model in a two-dimensional lattice 

icr i i i 



where n« = + n^, with = c ia Ci tCT . a(c) and a)(c^) are the phonon 
(electron) annihilation and creation operators. For simplicity we consider only 
the nearest-neighbor hopping tij = —t for i and j nearest-neighbor lattice 
sites. The above Hamiltonian defines the starting point of our flow. Since under 
the flow new interactions will be generated, in order to describe the evolution 
of the Hubbard-Holstein parameters, we have to write the Hamiltonian in a 
more general form 



H=T+V+G 

T = Y,u q : a\a q : + ^(e fe - /i) : cfcj : 



V 



k,T 



X! Vk,k',q '■ C k+q C k>-q C k' Tc k '■ 



2N 



k,k' ,q,T 



G = (9k,q a-q + 9*k+q,-q « 9 ) = 4+ ? C fe : 

V k,q,T 



(2) 
(3) 

(4) 

(5) 
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Figure 1 shows the diagrammatic representation of the el-ph interaction. In the 
following we assume the coupling g kq to be real, while u q and e k represent the 
phonon and electron dispersions respectively. V kj k',q is the e-e interaction. iV is 
the number of lattice sites. Colons indicate normal ordering of the operators 
with respect to the Fermi sea. The expression (2)represents the Hamiltonian 
at a generic value of the flow parameter I whereas Eq. 1 provides the initial 
values of the parameters at I — 0, i.e. gk, q (l = 0) = go, Vk,y,q{l — 0) = U, 
uj q {l = 0) = u , and e k = —2t(cos(k x ) + cos(k y )) (here we use a unit lattice 
spacing) . 




Fig. 1. Bare e-ph vertex. Streight lines represent electron propagators while wavy 
lines are for phonons. With our notation, in the vertex with an outgoing phonon 
the first index of g represents the incoming electron momentum, while in the vertex 
with incoming phonon it represents the outgoing electron momentum; in both cases 
the second index of g represents the opposite of the phonon momentum. 



2 Method 



The flow-equation method (jl2l : lira ) is based on an algebraic technique, the 
Jacobi method, to (block) diagonalize Hamiltonians. A series of infinitesimal 
rotations is performed such that at each step the non-diagonal elements get 
smaller while the others asymptotically approach a finite limit. If I is the pa- 
rameter ruling the flow, this technique performs a unitary transformation U (I) 
such that H(l) = U{l)H{l = 0)W(l). One can absorb all the evolution in the 
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parameters of the Hamiltonian generating a new H(l = oo). This transformed 
Hamiltonian has the property that, when expressed in terms of the original 
operators, it has the desired diagonal (or block-diagonal) form with renormal- 
ized couplings. The spectrum naturally stays unchanged. The evolution of the 
Hamiltonian is controlled by 



dJ ^ = WlH{l% (6) 



where the anti-hermitian generator rj(l) = = Proper identi- 

fication of the numerical coefficients in this operator equation yields the flow 
equations. In a finite Hilbert space, one can use a matrix representation of 
the H and r\ operators bringing Eq. 6 to a closed system of equations. On the 
other hand, working in second quantization, in general the commutator in Eq. 
6 will generate iV-body operators, which are not necessarily included in the 
starting Hamiltonian. In this case a truncation is needed at some N in order 
to obtain a closed form. For operators of conserved particles, this truncation 
is naturally justified in the low-density limit. On the other hand, for phonons, 
this truncation is only allowed within a weak-coupling approximation. 

The first step is to find a generator that produces the convergence of the flow 
to the desired form of the Hamiltonian. We split the Hamiltonian in H plus a 
part that has to be eliminated V. Wegner has shown that the following choice 
produces the desired flow 

BET 



17(0 = [H (l), V{1)] = [H (l), H(l)} = [H(l), V(l)}. (7) 



This method has been used to eliminate one kind of interaction (electron- 
phonon) in favor of an other effective (electron-electron) interaction (fl7h . The 
method was also adopted to select a specific channel in the electron-electron 
Hubbard interaction to see how this is dressed by all other channels in order 
to study the different possible instabilities of the model! 14 : Hit Il6h . Here, as 



a first step starting from the model Hamiltonian (2), in the weak correlation 
limit, we use this technique to eliminate the e-e interaction and study how 
the resulting effective e-ph interaction behaves. Since this use is somehow 
unconventional we will first test if the method is suitable to this purpose in a 
toy model with two sites and a truncated phonon Fock space and then proceed 
to the general case. Subsequently we will consider the (more standard fl?])) 
opposite procedure to eliminate the e-ph interaction in favor of an effective 
e-e interaction. 
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3 Effective Electron-Phonon interaction 



3.1 Two-site lattice 



As mentioned above, our first aim is to eliminate the e-e interaction term 
dressing the e-ph one. In order to achieve this result we use the generator 

r ] (l) = [T + G,V]. (8) 



In this case we are not properly trying to diagonalize the Hamiltonian but 
we want to eliminate some off-diagonal elements corresponding to the e-e 
interaction to see how other off-diagonal elements, the ones corresponding to 
the e-ph interaction, renormalize. Thus it is necessary to check if the generator 
given above yields a solution with vanishing V^k^q and finite effective coupling 
<7fc j9 . An analytical proof is beyond the purposes of this paper and therefore we 
will proceed constructively by analyzing a two-site Hubbard-Holstein model 
as a first non-trivial playground to apply this technique. 

One first remark is in order: The in-phase ionic displacement corresponds to a 
rigid shift of the two-site "lattice" and decouples from the electronic degrees of 
freedom. Therefore in momentum space the phonons at zero momentum can be 
disregarded and we only will consider phonons at momentum n. Accordingly 
only the e-ph couplings g^ q with q = n total transferred momentum will be 
considered. Specifically we consider a toy model model consisting in a two-site 
lattice with two electrons which can occupy a truncated Hilbert space made 
by the following S z — states: 

|1> = |ctT c ti| ) , |2> = \<$<#\0) , |3> = |44 T aO> , |4> = i4c^4 T |o> • 

The crucial simplification of this model is that the phonon Hilbert space is 
truncated to the subspace with only zero or one phonon. Of course this is 
a severe limitation if one aims to apply the model to the case of strong e-ph 
couplings, but is acceptable, instead, for the general analysis of the method and 
for the investigation of the weak e-ph coupling regime. We use the following 
single-electron dispersion 

e q = -2tcos(g). (9) 
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In matricial representation (2) takes the form: 



H(l) 



MO V 0fiA l ) 9n,A l ) -9irA l ) 



\ 



v ,oA l ) e o(0 

(0 9oA l ) 



9oA l ) ~9oA l ) 

a(0 -iw(0 



\-9nA l ) -9oA l ) - v o,nA l ) A (0 J 



(10) 



where 



e ff (Z = 0) = 2e 7T + U 
e o (Z = 0) = 2e + £7 
A(/ = 0) = e^ + e + f7 + w . 



(11) 
(12) 

(13) 



Since we aim to eliminate the e-e interaction and see how the e-ph interaction 
changes, we must implement a transformation to eliminate the matrix elements 
(1,2), (2,1), (3,4) and (4,3) of H. With the generator rj(l) given by Eq. 8 the 
flow equations read 



dl 

de (l) 
dl 



= -^9oA 1 )9ttA 1 ) V o,oA 1 ) + ^97tA 1 ) 2v o,7tA 1 ) 
+ 2V 0fi Al) 2 (eAl)-eo(l)) 

= -4:g Al)9nA l ) V o,oA l ) + ^A l f V ^A l ) 



+ 2V , Al) 2 (eo(l)-e 7r (l)) 

dX{l) =Ag Ai)9nAW,oAi) -2^(0 2 W(0 -w , v Ai)9«Ai? 

= -(e (Z) - e n (l)) 2 V ,oAl) - 2g Al) 2 Vo,oAl) 

+ 4:9oA l )9nA l ) V 0,7rA l ) ~ ^A l f V QflA l ) 

= -2g Al) 2 V^A l ) + 4goA l )9nA l )Vo,oAl) - ^A l f v ^A l ) 



dl 

dVofiAJ) 

dl 

dVo,irA l ) 

dl 
dl 

dgoAJ) 

dl 



(14) 

(15) 
(16) 

(17) 
(18) 



^(0^(0(2^(0 - e (Z) - A(Z)) + g ^(l)V^A^-Ul) + A(Z)) 



2#0,7T (Z)V ,0,7r(Z)V r o, 7r , 7r (Z) + fi' 7 r,7r(Z)K),0,7r(Z) 2 + 9nA l ) V 0,^A l f 



(19) 



: + A(0) + ^(0 W(0(-^(0 + 2e (Z) - A(Z)) 



2^(0^0,0,^(0^0,^^(0 + ^0,^(0^0,0,^(0 2 + ^(0^W(0 2 



(20) 



These equations can be solved numerically. As a consistency check we have 
compared the eigenvalues of both H(l = 0) and H(l = oo) for different values 
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of the starting couplings and obtained an exact match within our numerical 
precision. 



In the many-site case we cannot deal with the full form of the flow equations. 
Therefore it is useful to test which approximation can be done without gen- 
erating neither too big errors in the asymptotic values of the effective e-ph 
coupling nor a divergence in the case of degenerate levels. The simplest choice 
is to make a weak-U approximation stopping at first order in U. Aiming to 
the effective e-ph interaction, we can neglect the flow equations for the ener- 
gies because they would give higher order corrections when introduced in the 
flow-equations for the g's: 



de^l) _ de (l) _ d\(l) _ Q ^ 



dl dl dl 



dl 

+ 4<? o , 7 r(0<M0 W(0 - 2 <M0 2 tW0 (22) 

dV ^ d f ] = -2^(0 2 W(0 + 4<?o, 7r (0<M0^o,o, 7 r(0 

-2g^(l) 2 V ^(l) (23) 

dg^A 1 ) 



= <7o,,r(0tW(0(2e*(0 " eo(0 - A(0) 

+ 5^(0^0,^(0 (-^(0 + A(0) 

= 5o,.(0W(0(-eo(0 + A(0) 



dl 

+ ff 7 r, 7 r(0W(0(-e 7 r(0+A(0) (24) 

dgoA 1 ) 



dl 

+ gnAl)Vo,oAl)(-Z*(l) + 2e {l) - A(0) 



(26) 



We analyzed the behavior of the approximated effective e-ph couplings com- 
paring them to the exact case for different values of go and U G [0, 3], ujq = 
0.02. For the sake of convenience, we have used energy units such that t — 1/2. 
As far as ujq is concerned, we consider the most frequent adiabatic regime with 
a typical value corresponding to a hundredth of the maximum electron band- 
width. 



In figure (2) we show a comparison between approximated and exact asymp- 
totic values of gk, q for various values of U. In this example strong renormal- 
izations of g^ A at various U are present, both in the exact and in the approx- 
imated expressions. The approximated solutions are in good agreement with 
the exact ones only for values of electronic repulsion U smaller than t. Both 
these phenomena, the strong renormalizations of g and the sizable differences 
between the exact and the approximate results, are due to the big energy dif- 
ferences between the states with and without phonons, basically controlled by 
the energy-level difference of order t and to the small value of g we have been 
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Fig. 2. Exact (thick) and approximate (thin) numerical results for go j7T (green) and 
^^(red) vs U. Here we assume gk, q {l = 0) = 0.3, Uq = 0.02. 

using. As far as the second point is concerned, the approximated and the exact 
values are different because the vanishing of Vo,o,7r is controlled by t, whereas 
#o controls the vanishing of Vo i7rj7r , which for small values of go goes to zero 
slowly in the flow equations. The presence of g% is necessary in the equation 
for Vo.-tt.tt to get an asymptotically vanishing solution. Otherwise, due to the 
degeneracy of the levels A, Vo j7ri7r would not flow at all. The quadratic terms 
in U, which are missing in our approximated expressions, modify strongly in 
this discrete case the asymptotic values of g when U is not very small. 

To further clarify the physical interpretation of the various terms appearing 
in the exact asymptotic form of the Hamiltonian, we rewrite in real space 
the e-ph interaction. By disregarding the zero-momentum phonons, the e-ph 
coupling for a fixed value of the spin takes the form 



{g s Ai + g s A)6 + (gjA - gjA^)U 



(27) 



where 5 = c\c\ — c\c2, J = i(c\c2 — cjjci), A = a\ — a2 and 

90,tt + 9n,n 



95 



9.1 



(2* 



We explicitely checked that at very large U, gj grows monotonically with U 
while gs decreases vanishing asymptotically in U. This result naturally arises 
from the correlation-induced suppression of charge fluctuations. Intuitively, as 
U increases the charge fluctuations become stiffer and this reflects in a decrease 
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of the e-ph coupling that couples with the charge inbalance. Previous analyses 
already found that at small transferred momenta, the vertex corrections to g 
arising from the e-e interactions suppress the e — ph coupling by a factor (1 + 
Fq), where Fq is the standard Landau parameter entering the compressibility 
and is an increasing function of U (0; 0). Indeed the out-of-phase phonon A 
couples with the charge unbalance between the two sites. Since charge-density 
fluctuations are hindered by the local repulsion U, this is mirrored at large 
U in the suppression of the related coupling gg. On the other hand, the flow 
equations implement a transformation, which obviously keeps the spectrum 
unchanged. Since the energy scale U is present in the starting Hamiltonian, 
the same scale must appear in the transformed Hamiltonian as well, gg vanishes 
for large U's, and the interacting system keeps the scale U via the dressed e- 
ph interaction gj. In particular, in the large-?/ limit (which can be dealth 
with exactly in the t = limit) one can show that the local repulsion U 
is replaced by a nearest-neighbor attraction V ~ gj/^o- This leads to the 
relation gj ~ ^ujqU, which we numerically tested in the large-U limit. The 
appearance of a non-local e-ph coupling is also obtained by a standard Lang- 
Firsov trasformation 

c 9 
exp olf = ex P 

which is usually implemented to eliminate the linear local e-ph coupling in 
favor of an effective local e-e attraction —2g 2 /ujQ and a coupling between the 
phonons and the hopping term. By expanding to first order in the phonon 
fields this last term one can easily check that a linear non-local e-ph coupling 

gj = - tC - J2 ( c laCja - c] a c ia ) (at - a) - (a { - a,)) . 



is obtained. Then the specific value g = JUujq/2 also leads to the elimination 
of the e-e repulsion in the transformed Hamiltonian. Therefore for this very 
specific choice, the effect of the Lang-Firsof transformation is similar to the 
one generically obtained via our flow equations at large values of U, where 
the local e-ph coupling gg vanishes together with V, while gj(l = oo) tends to 

y/Uu /2. 1 



1 We also find in the Lang-Firsov-generated e-ph coupling a proportionality to the 
exponential factor exp[— (g/u>o) 2 ]. This factor is somewhat misleading since it should 
disappear in the effective four-fermion interaction as a result of the summation over 
intermediate multi-phononic states (see, e.g ., the effective superexchange coupling 
in the Hubbard-Holstein model of Ref. (jl9r ) or the effective Kondo coupling in the 
Anderson-Holstein model of Ref. H3)). Here the limitation of or 1 phonon states 
prevents this cancellation to occur. 



E 



QjA (Xi'. 



10 



Before going on analyzing the 2D continuum case, we reassert the two relevant 
indications of this section: The first is that it is indeed possible to eliminate 
the e-e interaction matrix elements obtaining a convergent solution for the 
effective e-ph couplings. The second indication is that the presence of an ef- 
fective phonon-mediated interaction introduces the quadratic terms in e-ph 
coupling in the flow equations of the e-e coupling, allowing to get a conver- 
gent asymptotic solution with a vanishing effective V already at first order in 
U. 

3.2 Effective e-ph interaction 

So far we discussed the generalities of the method and showed the possibility 
to have a convergent solution while eliminating some matrix elements in favor 
of others, using the generator (8) within our toy model. Here we proceed to 
the calculation of the effective e-ph coupling in the more general case of a two 
dimensional continuum. In this work, for simplicity we limit our analysis to 
the first order in the bare e-e coupling so that the flow equations are obtained 
from the following commutators 

[[G,V],T] + [[T,V],G]+o(U 2 ) 
[[T,V],T] + [[G,V],G]+o(U 2 ) 
[[G,V],G]+o(U 2 ) 
[[G,V],G] + o(U 2 ). 

Since our aim is the determination of an effective e-ph coupling while elimi- 
nating U at first order, we do not consider the flow equations for the energies 
because they would add higher order terms in U when inserted in the equa- 
tions for the g's. Having neglected the U 2 terms, corresponds to the t 3> U 
and g 2 /t 3> U conditions. Thus the flow equations for gk, q (l) and V k ^ >q (l) read 

Ael k , tq V k , k/ , q (l) + «[[G,V],Gr (29) 

7 2 9p-q,q( l ) V k,P,q( l ) ( n p-q ~ n p) [^k,k',q + «fc,g] (30) 
P 

A€ k ,k',q = + € k'-q ~ ^k ~ Cfc 



dg k , q {l) 

dl 

dV k ,k',q(l) 
dl 

de k {l) 
dl 

duJgjl) 
dl 



dV k ,k',q{ 1 ) 

dl 
dg k , q {l) 

dl 

where 
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&k,q = Cfc+g — e fc + ^0 

and U [[G, V], G]" stands for "part of flow equations coming from the operatorial 
term [[G,V],G] n and makes it very difficult to get an analytical solution. 
Nevertheless, as in the toy model, this term is essential for those values of 
k = k, k' = k' and q = q such that Ae kk ,g = 0: For these momenta, without 
this term, Vj. & would not flow to zero, and Eq. 30 would be non convergent. 

The term [[G, V],G] is cumbersome, but for our purposes it will be enough 
to keep the simplest diagonal part and neglect additional non-diagonal ones 
with one of the three momentum indeces of V under summation 

"[[G,F],q"~-8 5o 2 ^J. (31) 



Working at first order in U, we have also replaced g{l) by its bare value g . 
The solution for 14^/^(7) is then 

IWO = U e-K^ +8 *°K (32) 



Again we can substitute g p - q , q in the sum in Eq. (30) with its initial value g 
and solve it 

/ u e -Wl, P , q +H\i _ i \ 

9k, q (l) = 9o I 1 - E +g ^2] ( n p-i - n p )[2Ae k ^ q + I (33) 



so that 



9k, q (oo) =g (l + U x(k, q, u )) . 



(34) 



where 



0(h n , i \ — 1 n v) [^Ae k .p,q + Ctp-q,q} 



To start with we assume e k + q — e k + ou = (elastic scattering) so that the 
k dependence in the asymptotic solution disappears and x(k,q,u>o) becomes 

x _ 1 ^ {n P -q - Up) {e p ^ q - e p - a; ) 
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If one further neglects the g$ in the denominator the asymptotic solution 
becomes 




where x(<L ^o) is the usual Lindhardt polarization bubble. It is worth noting 
that in the limit of small coo, since J2 P n ^" q ^ n J < 0, we have g g (oo) < go for 
every value of the transferred momentum q. On the other hand, the general 
expression for g allows also for positive correstions to g in the small-g region. 
Moreover in the small-cuo limit, since 




(here the density of states M is the usual static small-g limit of the Lind- 
hardt function) one obtains g = go(l — UN). This result coincides with the 
static limit of perturbation theory at first order in U, thereby supporting the 
reliability of our approximations. This result also perturbatively agrees with 
the finding that e — e interactions suppress the e — ph coupling by a factor 
(1+F S ) where at this order F S = MU is the standard Landau parameter 

entering the compressibility. 

Returning to the more general form of Eq. (34), in figure 3 we show the 
behavior of the relative change of the effective coupling g^ A . To simplify the 
numerical treatment we assume that electrons have a free particle dispersion 
(we will also assume ^- = 1): e q = q 2 , and the direct self energy terms 
are included in the chemical potential fi, which is taken to be unity. The 
constant density of states per spin is Af = where v is the volume of the 
unit cell. When one considers the case of one particle per site, v = 2ir. Notice 
that, consistently with the assumption of a free-electron band dispersion, no 
limit has been imposed to the bandwidth, thereby leading to a vanishing of 
the "polarization-like" functions x, x an d X m t ne limit of large transferred 
momenta q — > oo. 

The effective coupling is unchanged at q — 0. This is due to the difference 
( n p-q ~ n p) which vanishes in the numerator and is not compensated by the 
denominator (as it would instead happen in the Lindhardt function) because of 
the <7q factor, which makes the reduction of the e-ph coupling small at small g's. 
Conversely, upon increasing q the reduction of g becomes more pronounced. 
This occurs up to (vpq + u Q ) 2 ~ g%, when the x is we U approximated by x- 
This latter in two dimensions, (if lattice effects are neglected), is constant up 
to q — 2k F (and, with our units, x — = —1/2) and then vanishes. As a 
result g reaches a plateau and finally recovers its bare value at large g's. 
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Fig. 3. The relative variation 9k ' 9 g^\j 90 °f the e-ph coupling. Here 
Ek,q = £fc+<j — £g = 0, loq = 0.02. We recall that with our units, M = 1/2 



The momentum dependence given in Fig. 3 shows the important feature of 
the e-ph interaction in the presence of correlation of being peaked at small 
momenta. The reduction of g at large transferred momenta has been already 
found in previous treatments of the Hubbard- Holstein model with slave bosons 



lOTl. Hubbard operators (9) and quantum Monte Carlo ([ill), where the 
momentum scale at whi ch q decreses is set by kp- Within a two-dimensional 
Monte Carlo approach ([ill), it was found that the decrease of g is largest 
at intermediate IPs, and becomes smaller at large U. Of course this non- 
monotonic result is not found with Eq. (34) here because of the weak coupling 
approximation. 



Having obtained an effective e-ph interaction, which is reduced by the Hubbard 
repulsion U and acquires a momentum dependence, it is of interest to study 
its effects on the electronic dispersion. One should, however, notice that the 
present first-order approximation limitates our possibility of obtaining sub- 
stantial quantitative changes in the electronic properties with respect to the 
standard Engelsberg-Schrieffer results in Ref. (J23J). We calculated the electron 
spectral density A(k,u) = —(l/ir)ImG(k,uj), where the electron Green's func- 
tion contains the self-energy corrections at second order in g. Fig. 4 displays the 
dispersions of the spectral-function peaks for various values of U at ujq = 0.02 
and go = 0.3. The overall behavior is obviously similar to that found in (23). 
Specifically the e-ph coupling mixes the electron with the phonon and bends 
the low-energy part of the dispersion on the phonon dispersion at E ~ c^o- On 
the other hand at very large energy the electron recovers its bare dispersion. 
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Fig. 4. Dispersion of the peaks in the electronic spectral function for different values 
of U, and go = 0.3, uj$ = 0.02. We report uj vs where u> is the solution of 
Cj — ek = Re[E(k,uj)] (filled symbols). We also report the position of the maximum 
of A(k,u) for uj > loq (unfilled symbols). 



It is clearly apparent that the e-ph interaction produces a break in the electron 
dispersion, which is naturally more pronounced for larger e-ph coupling (i.e. 
smaller U). The presence of weak interaction U introduces an additional in- 
teresting feature: Upon increasing U, the quantity A oc g 2 / (u> t) is suppressed 
thereby reducing the electronic mass m ~ m (l + A) near the Fermi level. This 
is made visible in Fig.4 by the increased velocity of the electrons at low energy. 
This decrease of the mass upon increasing the interaction U is a generic fea- 
ture, which was also found within a DMFT treatment of the Hubbard- Holstein 
model (|2lh. We also notice that a difference in the low-energy and high-energy 
dispersions as well as the s-shaped dispersion of Fig. 5, (which displaces upon 
tuning the e-e interaction), are reminescent of the kinks (|25l 0) and of the 
breaks (Jirjh in the electronic dispersion of the cuprates experimentally ob- 
served along the (1,±1) directions and around the (±1, 0), (0, ±1) regions of 
the Brillouin zone respectively. These breaks are a rather generic feature of 
electrons coupled to collective modes (27; 2^; 29l: liT 



4 Effective electron-electron interaction 



An effective e-e coupling in the static limit can be obtained from the effective 
e-ph coupling of Eq. (34) considering the following two contributions diagram- 
matically depicted in Fig. 5 
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Fig. 5. e-e vertex from renormalized e-ph interaction. 

rreff 9kJJ = °°)9k'-q,q(l = OO) 9k+q,-q(l = Op) 9 k',-q(l = gg) / qq n 
V k,k',q ~ "J "J l dy J 



If one restricts the external momenta k, k' to the domain I such that the 
condition of energy conservation ty — e^-q = £k+ q ~^k is satisfied, one obtains 
a simple expression for V^jJ 



Veff _ 2 90 
V k,k',q ~ 

*9qU j2 ( n p-i ~ n p )[2(e k+q + e p _ 9 - e p - e k ) + (e p - e p _ g )] ^ 
fc, fc' G / 

In this way we obtain the first correction of order U/t to the purely phononic 
limit U e ff w —2qI/uoq 2 We also notice that, imposing the elastic-scattering 



2 The use of the perturbative diagrammatic approach in Fig. 5 implies the smallness 
of the e-ph coupling g <C t, since otherwise higher order phononic corrections (self- 
energy and vertex) should be considered. It is worth noticing that this limitation is 
instead not required in Section 4.2, where the elimination of the coupling V with 
the neglect of some commutators imposes conditions on the smallness of V coupling 
only. 
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condition e k >-. q — e k > = e k+q — e k = one retrieves the perturbative result 

V q ef J elas u c = -(2s 2 M)) (1 + Ux(q,u ) + U X (q, -u )) = -&9oM i 1 + 2 ^x(9, 0)) 
where \ is the modified "polarization" bubble of Eq. (36). 

We have obtained a reduced attractive interaction modified according to the 
incipient effect of correlation on the e-ph interaction. In order to get further in- 
sight on the interplay between electronic correlations and the e-ph interaction 
we turn to the study of the effective e-e coupling by reversing the previous pro- 
cedure: We use the flow equation technique to eliminate the e-ph interaction in 
favor of the effective V(k, k', q). The generator of the unitary transformation 
is 

V (l) = [:T: + :V:,:G:} (41) 
leading to the flow equations 



An 

— ^[[:T:,:G :},:T :] + [[: T:,:G:],:V :] 

+ [[:V:,:G:],:T-] + [[:V:,:G:],:V-l (42) 
[[■ T :,: G :],: G :] + [[: V :,: G :],: G :] (43) 

Here we do not consider the flow equations for the energies under the condition 
that (t, uj ) ^> g . Moreover, when inserted in the equations for the effective 
Vk t k', q , they would add higher-order corrections in g . In the equation for G 
we neglect terms of order U 2 go (the last commutator in the r.h.s.). In the 
equation for V we disregard terms of order UgQ (the last commutator in the 
r.h.s.). Both these approximations are valid under the condition t^> U . 

When rewritten in terms of couplings and with the above approximations, the 
flow equations read 



d 9k, q (l) 2 „ /,X 

^ = ~ a k,q 9k,q{ L ) 

- -^Y,9p-q,q{ l ) V Kp^ l )( 2a p-q,q + Ae k,p,q) {^p-q ~ Tip) (44) 
i V p 

dV k ,k>, q {l) , , \ n\ n\ 
— = - {a k '-q, q + a ktq ) g k >- q , q ( l ) 9k,q[l) 

- (a k +q,-q + OCk'-q) gk+q,-q(l)gk',-q(l) (45) 

where, again, a k>q = e k+q - e k + u> , Ae kjk > jq = e k+q + e k >- q - e k > - e k , and 
e k = —2t[cos(k x ) + cos(fey)]. The U = case was previously considered in 
Ref. (fl7h. In this limit, Eq. (44) and (45) can be solved without any further 
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approximation yielding 

9k,q(l) = 9o Exp{-a 2 kq 1} 



(46) 



for the e-ph coupling and 

t t /, \ 2 - a k'-q,q - ®k,q 2 ~ a k+q,-q ~ <*k' -q 

V k , k >, q (l = oo)= g — 2 — — + g — 2 —2 (47) 

a k'-q,q a k,q a k+q,-q a k',~q 



as the new effective e-e coupling. As in Ref. (17) Eq. (47) in the BCS channel^' 
—k) gives 



fc '" Ml j " 9o 2(e k -e k+q y + 2u* 9o 2(e k - e k+q )* + 2u; 2 



(e* - e fc +<?) 2 + 



(4* 



while in the case of scattering with energy conservation (e k + q — = —tk'-q + 
e fc /) it gives 

t r n \ 9o 9o 2 2o; , , 

Velastic{l = OO) = = O? j (49) 



As discussed in Ref. fll7[ ). the effective interaction in the BCS channel is always 



attractive, with a noticeable difference with respect to the resnlt of Erolich fi) 

In the generic U > case the equations (44) and (45) are not analitically 
tractable and one should treat them numerically. Here instead we aim to get 
an insight on the enegy scales involved in the phonons dressing the effective e- 
e interaction and therefore we consider an approximate treatment. We expect 
that the dependence of the effective couplings on the momenta of the fermionic 
legs is less important than the dependence on the transferred momentum. This 
is particularly reasonable if one is mostly interested in the physically relevant 
region of scattering processes around the Fermi surface. According to this 
expectation, we assume here that g^q mostly depends on the transferred mo- 
menta q and we neglect its k dependence. Under this assumption the electronic 
parts of the a's in Eq.(45) cancel leaving lj as the only scale and the k and k' 
dependence of V drops. With the further assumption that in Eq. (44) V can 
be replaced by its bare value U, Eqs. (44) and (45) become 



= -K q )9q(l) ~ ^ Efe - e P - q )(n P - q ~ n p ) (50) 
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^f = -4^(0 (51) 

In deriving the equation for g we exploited the obvious property Z) p (2u;o + 
(ffe — £k-q)(np-q — n p) — 0. Moreover to eliminate the k dependence from the 
r.h.s., we consider the average quantity a 2 , ~ (a% q ) over the momenta k 
around the Fermi surface. The resulting expression for (a 2 ) can be cast in the 
form 

(a 2 ) =t 2 f3 q + LU t lq + iul (52) 



where the functions f3 q and ^ q depend on the specific average procedure (e.g., 
from the width of the shell around the Fermi surface over which the average 
is carried out). In particular, the average could be taken strictly on the Fermi 
surface, thereby getting explicit expressions. In general one can notice that 
both (3 q and 7 g vanish for q — > and that (a^ ) is obviously a non-negative 
quantity. 

We have not reported here the cumbersome expressions of the 0(U 2 go) and 
0(g 2 U) terms, but we explicitely checked that under the conditions V k ^ q ~ U, 
adopted in the evolution for g and gk, q (l) ~ <7g(0> these additional terms vanish 
and no further corrections should be considered in the flow equations. There- 
fore within this approximation, the limitation of small U can be somewhat 
relaxed. 

The solution of Eq. (50) is 

g q (l) = g Exp{ [ - (t 2 [3 q + tuj^ q + u%) - —— ^(e p - e p _ q )(n p - q - n p )]/}(53) 

p 



Where we use the notation e p = 2te p . Here we note that the exponential 
argument is manifestly negative and it implies the vanishing of g q (l) for / — > oo. 
We can now calculate the effective electron-electron coupling 

V k , k , q (l = oo) = U - —— 2 J t f"°l _ V( r.(54) 

t 2 (3 q + tuo'jq + u$ + — E P (e P - tp- q ){n P - q - n p ) 



The appearance of lo only in the numerator is a consequence of the above- 
mentioned cancellation of the electronic scale in the a's strictly valid under 
the assumption that g q does not depend on k. The effective coupling has 
four non-negative contributions in the denominator, with purely phononic, 
purely electronic and mixed character: tu 2 , t 2 (3 q , tcoo'jq and 2Utr) q = ^ Y,p(tp — 
£p-q)( n p-q — n p) The last three expressions vanish at q — (0, 0) and one finds in 
this limit —2g 2 /uj as a correction to U. The presence of correlation generically 
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reduces this attractive contribution to Vk,k',q when q is finite. Depending on 
the relative importance of the various terms appearing in the denominator of 
Eq. 54 as determined by the parameters t, U and Uq, the effective interaction 
acquires different limiting forms. First of all we notice that the correct anti- 
adiabatic limit (uo/t), (u /U) — > oo is recovered 

Vwa(l = oo) = U-^. (55) 



which was also obtained in the limit q — > 0. When loq ^> t, a further momentum 
dependent correction •y q t/u is obtained. 

In the opposite limit, (ujq — > 0, but with A = 1g\/{u)§ 1) finite) we find that 
the whole contribution of the e-ph attraction vanishes for ujq — ■> 0. This is 
reasonable, since an instantaneous effective interaction cannot be affected by 



a static lattice deformation. We notice that in Ref. (|21l). where the effective U 
was also calculated within the DMFT approach, a phonon correction linearly 
vanishing with u> was also found. Numerically the correction involved the 
scales U ~ t and ujq. An effective U was also obtained by means of a slave- 
boson-Lang-Firsov treatment of the Hubbard-Holstein model (J22I ). where in 
the adiabatic limit t ^> Uq, it was found U e ff ~ U — g 2 /At for a square lattice 
in two dimensions. In our case, in the same adiabatic limit, our Eq. (54) results 
in 

Ueff^U-—^- ^. (56) 

ff t/3 q + 2U Vq t K } 



Our present crude approximate treatment therefore reproduces the results of 
the DMFT and of the Lang-Firsov-Slave-Boson approaches with an additional 
uo/t factor. This factor in this limit provides a substantial reduction of the 
phononic corrections at finite g's, while this reduction is less and less severe at 
small transferred momenta, owing to the vanishing of (3 q and r] q in this limit. 
On the other hand, going back to Eq. (51), one should remember that relaxing 
the assumption of k independence of g, an additional electronic scale coming 
from the a>k, q would appear together with loq. As long as Uq is larger than 
t, these corrections are unessential and the present approximation works well 
in the antiadiabatic limit. On the contrary, in the adiabatic limit and finite 
q these electronic corrections should appear in the numerator of Eq. (54). 
As a consequence U would be corrected by a term of order g 2 /t without any 
additional uio/t factor. This analysis therefore suggests that the flow equations 
correctly capture the proper energy scales and provide non-trivial momentum 
dependencies. 
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5 Conclusions 



This paper considered two main issues related to the interplay between e-e cor- 
relations and e-ph coupling. We first investigated the effect of the e-e Hubbard 
repulsion on the Holstein e-ph coupling, showing that this latter is generically 
suppressed. This finding is consistent with previous results based on general 
Fermi-liquid arguments (0), on large-N expansions based on the slave-boson 
(0; 0) or Hubbard-operator (0) techniques and on numerical quantum Monte 
Carlo approaches (|llh Remarkably, although we limited our discussion to the 
weak-correlation limit, we also find an agreement as far as the momentum 
dependencies are concerned. Specifically we find that the Hubbard repulsion 
more severely suppresses the e-ph coupling when large momenta are trans- 
ferred. 



As already mentioned in the introduction this finding could explain the lack 
of phonon features in transport properties (which mostly involve large trans- 
ferred momenta) while charge instabilities and (i-wave superconductivity can 
arise from (mostly) forward scattering and still be induced by a less suppressed 
e-ph coupling. The suppression of the e-ph coupling by e-e correlations is also 
found within the Dynamical Mean-Field Theory technique, which provides an 
exact treatment of the dynamics. However, although this technique can be 
applied for any strength of the e-e and e-ph couplings, it neglects the space 
correlations and provides poor informations on the momentum structure of 
the renormalized couplings. In this sense our work based on the flow-equation 
technique explores the correlation-induced renormalization of the e-ph cou- 
pling from a complementary point of view. 

A second outcome of our research is provided by the elimination of the e- 
ph coupling to obtain an effective e-e interaction. Remarkably our procedure, 
even though it has been discussed within crude approximations, produces a 
reliable effective model for correlated electron systems in the presence of the 
e-ph interaction. All the energy scales, namely uo,go,t and U are explicitely 
included and provide proper momentum dependence corrections in the differ- 
ent physical regimes in agreement with previous analyses carried out with the 
DMFT and with the lang-Firsov-Slave-Boson approaches. 

The phonons naturally introduce an effective attractive interaction, which re- 
duces the repulsive U. Again we find that the phonon-mediated interaction is 
more sizable at small transferred momenta. This allows to infer some prelimi- 
nary conclusions on the effects of phonons on the instabilities of the Hubbard 
model. This issue had previously been addressed within the flow-equation 
technique where a phase diagram was proposed for the various insta- 

bility lines as a function of U and doping x. Here we qualitatively remark that 
the phonons induce an additional attraction, such that the effective repulsion 
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U(q) is smaller at smaller momenta. Therefore, the instabilities occurring at 
large momenta (like antiferromagnetism or flux-phase) are only weakly af- 
fected by the phonon attractive part, while the instabilities occurring at small 
momenta (phase separation, Pomeranchuk, superconductivity) will occur at 
larger values of the bare repulsion parameter U. 
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